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ABSTRACT 

The amplitude of the gravitational radiation from an accreting neutron star 
undergoing polar magnetic burial is calculated. During accretion, the magnetic 
field of a neutron star is compressed into a narrow belt at the magnetic equa¬ 
tor by material spreading equatorward from the polar cap. In turn, the com¬ 
pressed field confines the accreted material in a polar mountain which is mis¬ 
aligned with the rotation axis in general, producing gravitational waves. The 
equilibrium hydromagnetic structure of the polar mountain, and its associated 
mass quadrupole moment, are computed as functions of the accreted mass, 
M a , by solving a Grad-Shafranov boundary value problem. The orientation- 
and polarization-averaged gravitational wave strain at Earth is found to be 
h c = 6 x 10 _24 (M a /M c )(l + M a 6 2 /8M c ) _1 (//0.6kHz) 2 (d/lkpc)” 1 , where / is the 
wave frequency, d is the distance to the source, b is the ratio of the hemispheric to 
polar magnetic flux, and the cut-off mass M c ~ 1O~ 5 M 0 is a function of the natal 
magnetic field, temperature, and electrical conductivity of the crust. This value 
of h c exceeds previous estimates that failed to treat equatorward spreading and 
flux freezing self-consistently. It is concluded that an accreting millisecond pul¬ 
sar emits a persistent, sinusoidal gravitational wave signal at levels detectable, in 
principle, by long baseline interferometers after phase-coherent integration, pro¬ 
vided that the polar mountain is hydromagnetically stable. Magnetic burial also 
reduces the magnetic dipole moment /j monotonically as /i oc (1 + 3M a /4M c ) -1 , 
implying a novel, observationally testable scaling h c (n). The implications for the 
rotational evolution of (accreting) X-ray and (isolated) radio millisecond pulsars 
are explored. 

Subject headings: accretion, accretion disks — gravitation — stars: magnetic 
fields — stars: neutron — X-rays: stars 


1. Introduction 


The current generations of resonant bar antennas and long baseline interferometers are 
capable of detecting gravitational wave signals from incoherent sources with strain ampli¬ 
tudes h exceeding ~ 10~ 20 and ~ 10~ 21 respectively at frequencies near / « 0.6 kHz (Schutz 
1999). The sensitivity can be improved by a factor (/r) 1 / 2 for periodic sources of known 
/ by integrating coherently for a total observing time r, as in hcirarchical Fourier searches 
(Dhurandhar et al. 1996; Brady et al. 1998; Brady & Creighton 2000). Several classes of 
sources with promising event rates have been identified in the kilohertz regime: coalescing 
neutron-star binaries (Phinney 1991), r-modes in young, hot neutron stars (Andersson 1998; 
Lindblom et al. 1998), and neutron stars in low-mass X-ray binaries whose crusts are de¬ 
formed by temperature gradients (Bildsten 1998; Ushomirsky et al. 2000). More generally, 
every isolated and accreting millisecond pulsar is potentially a kilohertz source, because the 
stellar magnetic field, which is usually not symmetric about the rotation axis, deforms the 
crust and interior hydromagnetically. It is presently believed that hydromagnetic deforma¬ 
tions are too small to produce gravitational waves detectable by the current generation of in¬ 
terferometers, with h rsj 10 31 (H/10 12 G) 2 (//0.6kHz) 2 (d/lkpc) 1 for an object with surface 
magnetic field H, situated at a distance d from Earth (Katz 1989; Bonazzola & Gourgoul- 
hon 1996). Magnetars, with B > 10 15 G, are invoked as a possible exception (Konno et al. 
2000; Ioka 2001; Palomba 2001), as are objects whose internal toroidal fields are many times 
greater than B (Cutler 2002), or whose macroscopically averaged Maxwell stress tensor is 
enhanced relative to uniform magnetization because the internal field is concentrated in flux 
tubes, e.g. in a type II superconductor (Jones 1975; Bonazzola & Gourgoulhon 1996). 

In this paper, we show that previous calculations materially underestimate the hydro- 
magnetic deformation of recycled pulsars, e.g. accreting millisecond pulsars such as SAX 
J1808.4—3658 (Wijnands & van der Klis 1998; Chakrabarty & Morgan 1998; Galloway et al. 
2002; Markwardt et al. 2002). During accretion, in a process termed magnetic burial , material 
spreads equatorward from the polar cap, compressing the magnetic field into a narrow belt 
at the magnetic equator and increasing the field strength locally while reducing the global 
dipole moment fi (Melatos & Phinney 2000, 2001; Payne & Melatos 2004), in accord with 
observations of low- and high-mass X-ray binaries and binary radio pulsars (Taarn & van de 
Heuvel 1986; van den Heuvcl & Bitzaraki 1995). In turn, the compressed equatorial magnetic 
field reacts back on the accreted material, confining it in a polar mountain which is mis¬ 
aligned with the rotation axis in general (Melatos & Phinney 2000, 2001). The gravitational 
cllipticity e (Brady et al. 1998) of the mountain, calculated here self-consistently, can ap¬ 
proach e ~ 10“', materially greater than for an undistorted dipole [e ~ 6x 10“ 13 (.B/10 12 G) 2 , 
e.g. Katz (1989)] due to the enhanced stress from the compressed field. A mountain of this 
size generates gravitational waves at levels detectable by the current generation of interfer- 
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ometers and can explain the observed clustering of the spin periods of accreting millisecond 
pulsars through the stalling effect discovered by Bildsten (1998). 

In §2, we review the physics of magnetic burial and calculate the hydromagnetic struc¬ 
ture of the polar mountain as a function of accreted mass M a by solving an appropriate Grad- 
Shafranov boundary value problem, connecting the initial and final states self-consistently, 
and treating Ohmic diffusion semiquantitatively (Melatos & Phinney 2001; Payne & Melatos 
2004). In §3, we predict the orientation- and polarization-averaged gravitational wave strain 
h c (M a ) and compare it against the sensitivity of the Laser Interferometer Gravitational- 
Wave Observatory (LIGO). We also combine h c (M a ) with fi(M a ) from previous work (Payne 
& Melatos 2004) to deduce a novel, observationally testable scaling h c {p) for accreting mil¬ 
lisecond pulsars (Melatos & Phinney 2000). In §4, we discuss the implications of magnetic 
burial for the sign of the net torque acting on an X-ray millisecond pulsar and the evolution 
of h c and /i after accretion ends (e.g. in radio millisecond pulsars). Our conclusions rely on 
the assumption that the polar mountain is not disrupted by hydromagnetic (e.g. Parker) 
instabilities; the justification of this assumption is postponed to future work. Deformation 
by magnetic burial also occurs in white dwarfs (Katz 1989; Heyl 2000; Gumming 2002). 


2. Polar magnetic burial 

2.1. Equilibrium hydromagnetic structure of the polar mountain 

The equilibrium mass density p(x) and magnetic field intensity B(x) in the polar moun¬ 
tain are determined by the equation of hydromagnetic force balance, 

Vp + pV0 — (47 t) _ 1 (V x B) x B = 0 , (1) 

supplemented by (i) the condition V • B = 0, (ii) the gravitational potential 0 (with | V0| ~ 
GM^/Rl, where M* and R * are the stellar mass and radius, because the mountain is much 
shorter than R *), and (iii) an equation of state for the pressure p (isothermal for simplicity, 
i.e. p = CgP, where c s is the sound speed). We define spherical polar coordinates (r, 9, 0), such 
that 6 = 0 coincides with the magnetic axis before accretion, and we seek solutions to (1) that 
are symmetric about this axis, such that B can be constructed from a scalar flux function 
0(r, 9) according to B = (r sin d)^ 1 V0 x eLIpon resolving (1) into components parallel and 
perpendicular to B, we arrive respectively at the barometric formula p = c~ 2 F{pp) exp[— (0 — 
0o)/Cg], with 0 O = GM*/R *, and the Grad-Shafranov equation describing cross-field force 
balance (Hameury et al. 1983; Brown & Bildsten 1998; Litwin et ah 2001; Melatos & Phinney 
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2001; Payne & Melatos 2004), 

3 2 fi> sin 9 <9/1 3fi>\ 

dr 2 r 2 39 ysind 39 J 


—47rr 2 sin 2 9F'( / ip) exp[—(/> — 4>o)/cl ■ 


( 2 ) 


It is customary to choose an arbitrary (albeit physically plausible) functional form of 
Fffi) in order to solve (2) (Uchida & Low 1981; Hameury et al. 1983; Brown & Bild- 
sten 1998; Litwin et al. 2001). However, this approach leads to an inconsistency. In 
the perfectly conducting limit, material is frozen to magnetic field lines. Hence the mass 
dM = 2nd’ip J ds p[r(s), 0(s)]r sin 9\ V^| _1 enclosed between two adjacent flux surfaces if and 
ip + d'i/) (where the integral is along a field line) equals the mass enclosed before accretion plus 
any mass added during accretion, without any cross-field transport of material (Mouschovias 
1974; Melatos & Phinney 2001; Payne & Melatos 2004). To calculate if correctly, one must 
specify dM/dfi according to the global accretion physics, with F(if) following from 

F W = {/ds rsindlV^I-exph^-^/c, 2 ]} ; (3) 

otherwise, if F(ifi) is specified, dM/dfi> changes as a function of M a in a manner that in¬ 
consistently leads to cross-field transport. The functional form of F(fi>), determined self- 
consistently through (3), changes as M a increases. This property, newly recognized in the 
context of magnetic burial (Melatos & Phinney 2001; Payne & Melatos 2004), has an im¬ 
portant astrophysical consequence: it produces a greater hydromagnetic deformation than 
predicted by previous authors, because polar (accreting) and equatorial (nonaccreting) flux 
tubes maintain strictly separate identities through (3), without exchanging material, and 
hence the equatorial magnetic field is highly compressed. 

We solve (2) and (3) simultaneously for if(r,9) and p(r,9 ) subject to the line-tying 
(Dirichlet) boundary condition ip(R*,9) = sin 2 9 at the stellar surface, such that the 
footpoints of magnetic field lines are anchored to the heavy, highly conducting crust (Melatos 
& Phinney 2001; Payne & Melatos 2004). We adopt a mass-flux distribution, M(if) = 
\M a { 1 — eW/ba)(i _ ^ that embodies the essence of disk accretion, namely that 

the accreted mass is distributed rather evenly within the polar flux tube 0 < if < with 
minimal leakage onto equatorial flux surfaces where /y = n/2) denotes 

the hemispheric flux and is the flux surface touching the inner edge of the accretion disk 
(radius Rfi). Payne & Melatos (2004) verified that the results do not depend sensitively on 
the exact form of only on the ratio /; a /?/;*. 

The boundary value problem (2) and (3) with surface line-tying can be solved ana- 
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lytically, in the small-M a limit, by Green functions 1 after approximating the source term 
on the right-hand side of (2) to sever the coupling between (2) and (3) (Payne & Melatos 
2004). It can also be solved numerically by an iterative scheme that solves the Poisson 
equation (2) for a trial source term by successive overrelaxation then updates from 

(3) by integrating numerically along a set of contours (including closed and edge-interrupted 
loops) (Mouschovias 1974; Payne & Melatos 2004). The numerical scheme is valid for ar¬ 
bitrarily large M a , although convergence deteriorates as M a increases. Typically, we use a 
256 x 256 grid and 255 contours (linearly or logarithmically spaced) and target an accuracy 
of |At/?/?/?| < 1CT 2 after ~ 10 3 iterations. We rescale the r and 6 coordinates logarithmically 
in the regions where steep gradients develop, e.g. 9 ~ n/2 (Payne & Melatos 2004). 


2.2. Ohmic diffusion 

The structure of the polar mountain evolves quasistatically, over many Alfven times, in 
response to (i) accretion, which builds up the mountain against the confining stress of the 
compressed equatorial magnetic field, and (ii) Ohmic diffusion, which enables the mountain 
to relax equatorward as magnetic field lines slip through the resistive fluid. The competition 
between accretion and Ohmic diffusion has been studied in detail in the context of neutron 
stars (Brown & Bildsten 1998; Litwin et al. 2001; Gumming et al. 2001) and white dwarfs 
(Gumming 2002). In these papers, steady-state, one-dimensional profiles of the magnetic held 
are computed as functions of depth, from the ocean down to the outer crust, and the Ohmic 
(fd) and accretion (f a ) time-scales are compared, under the assumptions that the magnetic 
held is battened parallel to the surface by polar magnetic burial, the accreted material 
is unmagnetized, and the accretion is spherical; ‘the complex problem of the subsequent 
spreading of matter’ is not tackled (Gu mmin g 2002). The held penetrates the accreted layer 
if the accretion rate satishes M a < O.lMEdd, where Msdd is the Eddington rate, but it is 
screened diamagnetically (i.e. buried) if M a > O.lMEdd, such that the surface held is reduced 
~ (M a /0.002M E dd)-f°ld relative to the base of the crust (Gumming et al. 2001). 2 

In the regime fd 1$ f a , Ohmic diffusion outpaces accretion. As material is added, it 
does not compress the equatorial magnetic held further; instead, it diffuses across held gra- 


1 The Grad-Shafranov operator is not self-adjoint. 

2 SAX J1808.4—3658 is presumed to possess an ordered magnetic field at its surface because it pulsates 

(Wijnands & van der Klis 1998; Chakrabarty & Morgan 1998). Due to the low accretion rate, M a « 
lO - 11 M 0 yr -1 , either the Held has penetrated to the surface by Ohmic diffusion (Cunnning et al. 2001), or 
polar magnetic burial and equatorward spreading have not proceeded far enough (Payne & Melatos 2004). 
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clients and distributes itself uniformly over the stellar surface. Hence the polar mountain 
(i.e. the asymmetric component of p ) ‘stagnates’ at the structure attained when t d ~ t a . 
The accretion time-scale is defined as t a = M a /M a . The Ohmic diffusion time-scale is given 
by t d = 4ttctL 2 /c 2 , where a denotes the electrical conductivity and L = (\'ip\/\ V^Dmin is 
the characteristic length-scale of the steepest field gradients; L reduces to the hydrostatic 
scale-height in the one-dimensional geometry employed in earlier work (Brown & Bildsten 
1998; Cumming et al. 2001; Cumming 2002) but is dominated by latitudinal gradients here. 
Following Cumming et al. (2001), we assume that the electrical resistivity is dominated by 
electron-phonon scattering in the outer crust (M a > 1O~ 1O M 0 ), as expected if the crustal 
composition is primordial, although electron-impurity scattering may dominate if the prod¬ 
ucts of hydrogen/helium burning leach into the crust. In the relaxation time approxima¬ 
tion, with all Coulomb logarithms set to unity, the electron-phonon conductivity is given by 
a = 7.6 x 10 26 (p/10 11 gcm _3 )(i/ c /10 16 s _1 ) _1 (m e ff/m e ) _1 (H/2) _1 s _1 , where A and Z denote 
the mean molecular weight and atomic number per electron, m e g is the effective electron 
mass, u c = 1.2 x 10 18 (T/10 8 K) s^ 1 is the electron-phonon collision frequency, T is the tem¬ 
perature of the crust, and p = 6.2 x 10 10 HZ -1 (M a /10 -5 M o ) 3 / 4 gem 3 is the density at the 
base of the accreted layer (Brown & Bildsten 1998; Cumming et al. 2001). Ohmic diffusion 
therefore arrests the growth of the polar mountain for 


IW| 


< 4.2 x 10 2 Z 1/2 


Ma 

10- 5 Mq 


1/8 ( t \ l/2 fm eS y' 2 ( M a y 1/2 

U 8 kJ UeJ ly/EdJ cm 


( 4 ) 


corresponding to the condition t d < t a . The left-hand side of (4) is computed directly from 
the numerical solution, by scanning over the grid, or from the approximate analytic solution 
in §2.3; it depends implicitly on M a . We denote by M d the minimum accreted mass for which 
(4) is satisfied. Note that t d /t a is constant with depth for electron-phonon scattering but 
decreases with depth for electron-impurity scattering (Brown & Bildsten 1998; Cumming 
et al. 2001). 

Several second-order effects, neglected in (4), are postponed to future work. The Hall 
conductivity vanishes in plane-parallel geometry (Cumming et al. 2001) but more generally 
amounts to a fraction 1.8 x 10~ 3 (H/10 8 G)(z/ C /10 18 s” 1 ) _1 (mefj/m e ) _1 of cr (Cumming et al. 
2001, 2004), i.e. the electron cyclotron frequency divided by u c . It may therefore dominate 
near the equator, where the compressed field can reach B ~ 10 15 G. In addition, the process 
of thermomagnetic drift is negligible in the ocean and outer crust, but it can dominate in 
the thin hydrogen/helium skin overlying the ocean (Geppert & Urpin 1994; Cumming et al. 
2001). Neither the Hall nor the thermomagnetic drifts are diffusive; they do not smooth 
out field gradients as cr does. Indeed, the Hall drift tends to intensify field gradients by 
twisting the magnetic field in regions where the velocity of the electron fluid is sheared. The 
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time-scale of this process is sensitive to the field geometry (as well as the radial profiles of 
the electron density and clastic shear modulus in the crust) (Cumming et al. 2004); it will 
be interesting to evaluate it for the distorted held produced by magnetic burial. 


2.3. Mass quadrupole and magnetic dipole moments 

Figure la depicts the density profile of the polar mountain (dashed p contours) and the 
haring geometry of the magnetic held (solid ^ contours), for M a = 1O~ 5 M 0 and t/b/t/y = 
lO^ 1 . The accreted material spreads equatorward under its own weight until its advance 
is halted near the equator by the stress of the compressed (and hence amplified, by hux 
conservation) magnetic held. Figure 16 illustrates this balance of forces; the Lorentz force 
per unit volume, (47 t) _1 (V x B) x B (solid contours), peaks at the boundary of the polar 
hux tube that receives accreted material, while the local magnetic held intensity, |B| (dashed 
contours), is amplihed «4 x 10 3 times above its initial polar value. 

In Figure 2 a, we plot the gravitational cllipticity e = \I\ — h\/h, where I\ and I 3 
denote principal moments of inertia, as a function of M a for several values of The 

reduced mass quadrupole moment is then given by 2eI zz /3, with I zz « 0.4 M*R 2 (Shapiro 
& Teukolsky 1983; Bonazzola & Gourgoulhon 1996). There is good agreement between the 
numerical and analytic results for ^a/^* = 0.3, where the gradients are manageable and the 
code converges reliably. The analytic results follow from a Green function analysis in the 
small-M a limit, where the source term on the right-hand side of (2) is evaluated for a dipole, 


yielding (Payne & Melatos 2004) 

ip(x,y) =^*{ l -y 2 ){ l + x / a Y 1 [ l - (M a /M c )6 2 yg(y)f(x)\ , (5) 

p{x,y) = M a ab[2nR 3 J(b)]- l [y 2 + (1 -y 2 )(MjM c )b 2 yg(y)f(x)] 1/2 (6) 

x exp [—6(1 - y 2 ) - x + (. M a /M c )b 3 yg(y)f(x)\ , 

P = 1 + 3M a /4M c ) 1 , (7) 

e = (5M a /2M*)(l - 3/26)(1 + M a 6 2 , (8) 

to leading order in M a /M c , with 

M c = GM^l/AciRl . (9) 

We define the dimensionless quantities x = a(r/i2* — 1), y — cos 6, f(x) = 1 — exp(— x), 


g(y) = exp[—6(1 — y 2 )], a = GM*/c 2 R *, and 6 = Typically, the hydrostatic scale 

height is small compared to /?*, as is the polar fraction of the hemispheric flux, implying a 



1 and b 1 respectively. The critical accreted mass M c depends on the natal magnetic flux 
and crustal temperature through (9), with M c /M 0 = 3 x 1O _4 ('0*/1O 24 Gcm 2 ) 2 (T/10 8 K) -2 . 

The e(M a ) relation (8) associated with magnetic burial is new, while the /x(M a ) relation 
(7) is closely related to the phenomenological scaling invoked by Shibazaki et al. (1989) to 
model observations of binary and millisecond radio pulsars (Taam & van de Heuvel 1986; 
van den Heuvel & Bitzaraki 1995). It is important to note that magnetic burial affects e and 
/j in different ways. The accreted mass above which e saturates, 8M C 5~ 2 , and the saturation 
ellipticity, e max = 20(M C /M*)M 2 (1 — 3/2 b), are inversely proportional to b 2 (for b S> 1), 
whereas the accreted mass above which /i is screened, 4M c /3, is independent of b. This 
theoretical prediction conforms with observations in two key respects. First, most accreting 
millisecond pulsars have b > 10 (Litwin et al. 2001) and hence e max < 10 - ' from (8), consis¬ 
tent with the upper limit e < 10~ 5 inferred from spin down (Dhurandhar et al. 1996; Brady 
et al. 1998) and the failure of bar antennas and interferometers to detect gravitational waves 
so far (Schutz 1999). Contours of e max are plotted in the plane in Figure 2b. Second, 
the floor magnetic moment of recycled neutron stars, which is observed to be ‘universal’, is 
given by /x m ; n ~ 10 26 Gcm -3 (for M a ~ 1O _1 M 0 ) from (7). Theoretically, it is independent 
of b (and hence M a ) in the regime t d > t a , and weakly dependent on M a in the regime t d < t a 
(Shibazaki et al. 1989; Melatos & Phinney 2001). 

The dipole and quadrupole moments saturate at /t(M d ) and e(M d ) respectively when 
Ohmic diffusion dominates (M a > M d ; see §2.2). From (4), (5), and (9), we obtain 


Md 

M 0 


= 3.4 x 10" 7 Z- 4/9 




O.l^b 


16/9 


T 


10 8 K 


- 20/9 


i>* 


10 24 G cm 2 


16/9 


Mo 


4/9 


M, 


( 10 ) 


Edd 


Contours of M d in the M a -T plane are displayed in Figure 3. 

The distorted magnetic field in Figure la may be disrupted by interchange, Parker, and 
doubly diffusive instabilities wherever the local field strength exceeds B r-u 10 10 G (Gumming 
et al. 2001; Litwin et al. 2001; Payne & Melatos 2004). However, more work is required 
to settle this issue; existing stability calculations are linear and plane-parallel, unlike the 
situation in Figure la. The equilibrium field may not be disrupted completely if the insta¬ 
bility saturates promptly in the nonlinear regime; for example, the interchange instability is 
inhibited topologically by the line-tying boundary condition, which constrains the mobility 
of closely packed flux tubes. 3 In a separate effect, Payne & Melatos (2004) proved analyt- 


3 Indeed, the numerical solution in Figure la represents the endpoint of a convergent sequence of iterations 
in a relaxation scheme. Convergence is cited by Mouschovias (1974) as evidence of stability because the 
relaxation process ‘mimics’ (albeit imperfectly) the true time-dependent evolution. 
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ically that the magnetic field develops bubbles for M a > 1O _4 M 0 ; the source term on the 
right-hand side of (2) creates flux surfaces oc oc M a that are disconnected from the 

star (-0 < 0, V 7 > '0*)- Further work is required to determine how the bubbles evolve in the 
large-M a regime (e.g. M a ~ 1O _1 M 0 ), where the theory in §2.1 breaks down. 


3. Gravitational radiation 


3.1. Detectability 


The mass distribution resulting from polar magnetic burial is not rotationally symmetric 
in general; the principal axis of inertia e 3 (i.e. the dipole axis of the natal magnetic field) 
is inclined at an angle a to the rotation axis. Gravitational waves are emitted by the 
deformed star at the spin frequency Q = f/2 and its first harmonic 2 Q = f (unless a = 
7t/ 2, when there is no f/2 component), and the spin-down luminosity is proportional to 
(16 sin 2 ol Tcos 2 a) sin 2 a (Shapiro & Teukolsky 1983; Bonazzola & Gourgoulhon 1996). Upon 
averaging over a, polarization, and orientation (the position angle of the rotation axis on 
the sky cannot normally be measured), one can define the characteristic gravitational wave 
strain h c = (1287r 4 /15) 1//2 G'/ 22 / 2 e/(dc 4 ) (Brady et al. 1998), which reduces to 


h c = 7.7 x 10 


-19 


M a 

M* 


1 - 


2b 


1 + 


MJj 2 

8 Mr 


-1 


/ 


0.6 kHz 


d 


1 kpc 


-i 


(ii) 


upon substituting (8). Polar magnetic burial therefore generates gravitational radiation 
whose amplitude h c ~ 6 x 10~ 26 (for typical parameters M a > M c ~ 10 -5 M* and b = 30) 
is ~ 10 5 times greater than that produced by the natal, undistorted magnetic dipole, h c ~ 
10 -31 (.B/10 12 G) 2 (//0.6 kHz) 2 (d/l kpc) -1 (Katz 1989; Bonazzola & Gourgoulhon 1996), due 
to the enhanced Maxwell stress from the compressed equatorial magnetic field. The self- 
consistent form of F(^) defined by (3) is needed to calculate this stress properly; cf. Melatos 
& Phinney (2001). 

Figure 4a is a plot of h c versus / for 10 -8 < M a /M 0 < 10 -2 and 10 < b < 10 2 . 
The sensitivity curves for initial LIGO ( L) and advanced LIGO (AL) are superimposed, 
corresponding to the weakest source detectable with 99 per cent confidence in r = 10 7 s of 
integration time, if the frequency and phase of the signal at the detector are known in advance 
(Brady et al. 1998; Schutz 1999). The plot suggests that the prospects of detecting objects 
with M a > 1O -5 M 0 are encouraging, a point first made in the context of magnetic burial by 
Melatos & Phinney (2000). Such objects include the five accreting X-ray millisecond pulsars 
discovered at the time of writing, SAX J1808.4—3658, XTE J1751—305, XTE J0929—314, 
XTE J1807—294, and XTE J1814—338 (Wijnands & van der Klis 1998; Chakrabarty & 
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Morgan 1998; Galloway et al. 2002; Markwardt et al. 2002; Strohmayer et al. 2003; Campana 
et al. 2003). Detectability is facilitated if the time-dependent Doppler shift from the binary 
orbit is known well enough to be subtracted, but this is only practical for some objects of 
known /, not for an all-sky search (Brady et al. 1998). 

The characteristic gravitational wave strain h c is a lower limit in an important sense: it 
contains the assumption that, after averaging over polarization and orientation, the bulk of 
the gravitational wave signal is emitted at / = 2D. This is true when the observer views the 
system along its rotation axis and measures h + oc h x oc e 2?fM sin 2 a in the two polarizations. 
It is not true when the observer views the system perpendicular to its rotation axis and 
measures h + oc e 2 ' lQ,t sin 2 a and h x oc e lQt sin a cos a, for example (Shapiro & Teukolsky 1983; 
Bonazzola & Gourgoulhon 1996). One certainly expects, as a matter of chance, to observe 
some systems in the latter orientation, or close to it. For these sources, the average quantity 
h c substantially underestimates the true wave strain (which is dominated by h x in the above 
example), especially if a is small. 4 (The observed pulse modulation indices of accreting 
millisecond pulsars imply a range of a values.) The effect on detectability is even greater 
when one takes into account the shape of the interferometer’s noise curve; for example, LIGO 
is 1.5-2.5 times more sensitive at D ~ 0.3 kHz than at 2D ~ 0.6 kHz. 


3.2. h c versus /x 

Polar magnetic burial is not the only mechanism whereby accreting neutron stars with 
M a > 1O~ 5 M 0 can act as gravitational wave sources detectable by long baseline interferom¬ 
eters. Crustal deformation due to temperature gradients (§4) is one of several alternatives 
(Bildsten 1998; Ushomirsky et al. 2000). In principle, we can distinguish between these mech¬ 
anisms by eliminating M a from (7) and (8) to derive a unique — and testable — scaling 
/?. c (yu) for magnetic burial that relates observable quantities only, unlike (11), which features 
M a (usually inferred from evolutionary models) (van den Heuvel & Bitzaraki 1995). The 
scaling, 

h c oc (1 - 3/26) [/x/(x/>*-R* - /x) + 6 2 /6] _1 , (12) 

is graphed in Figure 4 b for 3 < b < 100. The vertical segments of all the curves represent 
the regime M a -C M c , where /x retains its natal value while the ellipticity grows as M a . The 
horizontal segments represent the regime M a M c , where e saturates while [i decreases, 
with the turn-off mass and e max scaling as b~ 2 . The proportionality (12) does depend on b 


4 We note in passing that h c also underestimates the true wave strain substantially for precessing radio 
pulsars, e.g. PSR B1828—11, whose wobble angle is known to be small (< 3°) (Link & Epstein 2001). 
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and neither of which can be measured, but the relevant scalings are moderately weak 
(e.g. b oc M a 2//7 V*' /7 ). Consequently, the overall trend in Figure 4 b may emerge statistically, 
once many gravitational wave sources have been detected, provided that the range of neu¬ 
tron star magnetic fields at birth is relatively narrow. Drawing upon population synthesis 
simulations, Hartman et ah (1997) inferred that ~ 90 per cent of radio pulsars are born with 
10 23 ' 5 Gcm 2 < yd < 10 24 ' 5 Gcm 2 , but the existence of anomalous X-ray pulsars implies a 
wider range of r/y in a subset of the neutron star population (Regimbau & de Freitas Pacheco 
2001 ). 


4. Discussion 

In this paper, we report on two new results concerning the gravitational radiation emit¬ 
ted by accreting neutron stars undergoing polar magnetic burial. First, upon calculating 
rigorously the hydromagnetic structure of the polar mountain, by solving a Grad-Shafranov 
boundary value problem with the correct flux freezing condition connecting the initial and 
final states, we find that the cllipticity of the star materially exceeds previous estimates, due 
to the enhanced Maxwell stress exerted by the compressed equatorial magnetic held, with 
e ~ 20M c /M*6 2 for M a > M c ~ 1O _5 M 0 , as given by (9). The associated gravitational wave 
strain at Earth, h c = 6 x 10~ 24 (M a /M c )(l + M a 6 2 /8M c )~ 1 (//0.6kHz) 2 (d/l kpc) -1 , averaged 
over polarization and orientation, is detectable in principle by the current generation of long 
baseline interferometers, e.g. LIGO. (The wave strain at //2 = in one polarization can 
exceed h c for many orientations and a values.) Second, the stellar magnetic moment /x 
decreases as the polar mountain grows and spreads equatorward, implying a distinctive, ob¬ 
servable scaling h c (/x), displayed in Figure 46, which can be used to test the magnetic burial 
hypothesis. Accreting neutron stars are easier to detect as gravitational wave sources than 
other kilohertz sources like coalescing neutron star binaries: they are persistent rather than 
transient, the waveform is approximately sinusoidal, and, if X-ray pulsations or thermonu¬ 
clear burst oscillations are observed in advance (e.g. SAX J1808.4—3658), extra sensitivity 
can be achieved by integrating coherently. 

The range of e predicted by the theory of magnetic burial is consistent with that in¬ 
voked by Bildsten (1998) to explain the clustering of spin frequencies (within 40 per cent of 
0.45 kHz) of weakly magnetized, accreting neutron stars (Chakrabarty et al. 2003) in terms 
of a stalling effect where the gravitational wave torque (oc fl 5 ) balances the accretion torque 
(oc M a R^ 2 oc Ma^'jU 2 / 7 ). Importantly, the theory of magnetic burial predicts that e increases 
monotonically with M a — a key precondition for the stalling effect to operate properly, oth¬ 
erwise the stall frequency G would not be a stable fixed point. In an alternative scenario, 
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Bildsten (1998) and Ushomirsky et al. (2000) attribute the requisite mass quadrupole mo¬ 
ment, 2 x 10”'/ 22 (M a /M Edd ) 1 / 2 (//0.6 kHz)” 5 / 2 , to lateral temperature gradients in the outer 
crust, which induce gradients in the electron capture rate and hence p. 

Magnetic burial, acting in concert with the stalling effect, influences the rotational 
evolution of an accreting neutron star in two observationally testable ways. First, during the 
late stages of accretion, the instantaneous net torque is zero, because the gravitational wave 
and accretion torques balance (Bildsten 1998), but the average torque (on the time-scale f a ) 
is effectively negative, because the gravitational wave torque increases with M a as the polar 
mountain grows; that is, the instantaneous stall frequency fl oc M a ^ 35 p,(M a ) 2 / 35 e(M a )” 2 / 5 
decreases with M a . This effect ought to be detectable by X-ray timing experiments in 
progress (Galloway et al. 2002), although it may be masked if M a fluctuates stochastically 
on the time-scale f a , as commonly happens. Second, magnetic burial predicts a distinctive 
evolutionary relation between p and fh During the early stages of accretion, before the 
star is spun up to the stall frequency, D increases while p decreases due to burial, with 
fl oc /?a / 2 M a oc M a / 1” 5 / 7 in the regime p oc M” 1 . However, during the late stages of 
accretion, ff and p both decrease as explained above. Hence there exists a maximum spin 
frequency Q/2n < 1 kHz bounding the population of accreting millisecond pulsars in the Q-p 
plane, and objects follow A-shaped evolutionary tracks in that plane. 

Once accretion ends, do we expect to detect the neutron star as a radio millisecond 
pulsar? There are two arguments against this from the perspective of magnetic burial. 
First, the gravitational wave spin-down time, t g ~ 5 x 10 7 (e/10” 7 )” 2 (ff/2krads” 1 )” 4 yr, is 
shorter than the observed age, so the neutron star rapidly brakes below the radio pulsar 
death line (O/27 t < 0.1 kHz) and is extinguished as a radio source. Note that this occurs 
no matter what generates the quadrupole moment inferred from the stalling effect. Second, 
this rapid spin-down is interrupted if the stellar deformation relaxes quickly (compared to 
t g ), for example if the buried, polar magnetic field is resurrected on the Ohmic time-scale f d 
of the outer crust, which satisfies f d ~ £a = 6 X 10 7 (M a /M 0 )(M a /M Edd ) 1 yr in the regime 
M a > 1O” 5 M 0 . However, the buried field is resurrected in stages; the steepest gradients are 
smoothed out (L oc t l J 2 ), but the natal held is not fully restored over the typical lifetime of 
a millisecond pulsar, because one has p oc L oc t]/ 2 and hence B ~ 10 9 ” 10 G after ~ 10 9 yr 
(Gumming et al. 2001; Mclatos & Phinney 2001). This scenario leaves the radio millisecond 
pulsars with the lowest fields unexplained unless they have low fields initially. 
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Fig. 1.— (a) Numerical calculation of the hydro magnetic structure of the accreted layer, 
showing contours of (solid) and p (dashed). (6) Contours of Lorentz force density (solid) 
and magnetic held intensity (dashed). Parameters: M a = 1O~ 5 M 0 , b = 10. The contours are 
at fractions r] = 0.8, 0.6, 0.4, 0.2, 10~ 2 , 10~ 3 , 10~ 4 , 10 -5 , 10 -6 , and 10 -12 of the respective 


maxima. 
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Fig. 2.— (a) Ellipticity e = 1 1\ — h\/h as a function of accreted mass M a , for 6 = ^*/^ a = 3, 
10, 30, 100. The curves are theoretical, based on the expressions (5)-(9). The points are 
numerical results for b = 3. Parameters: = 10 24 Gcm 2 , T = 10 8 K. (6) Contours of 

maximum ellipticity e max = 20M C /(M*6 2 ), from (7) and (9), as a function of total magnetic 
flux -0* and crustal temperature T. Contours are spaced logarithmically in the range 10" - 9 < 
M c /Mq < 10 -4 . In both panels, the regime td > t a is considered to highlight the effects of 
burial in isolation. 
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Fig. 3.— Contours of the minimum accreted mass for which Ohmic diffusion arrests 
the growth of the polar mountain (O < t a ). The plot shows the variation of Md with 
accretion rate M a (in units of MEdd) and crustal temperature T (in kelvins), in the range 
10~ 8 < M d /M Q < 10~ 5 , with i/j = 0.1, Z — 1, and = 10 24 Gcm 2 . M d is calculated 
from (4), (5), and (9). 
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Fig. 4.— (a) Polarization- and orientation-averaged gravitational wave strain h c as a func¬ 
tion of wave frequency /. The initial and advanced LIGO noise curves (solid) correspond 
to detection with 99 per cent confidence after 10's of coherent integration (Schutz 1999). 
Theoretical curves are shown for b = 10, M a /M 0 = 10^ 2 , 10 -4 , 10” 6 (dashed, top to bottom) 
and b = 10 2 , M a /M 0 = 10 -4 , 10 -8 (dotted, top to bottom). ( b ) Gravitational wave strain 
h c versus magnetic moment /U (scaled to the natal magnetic moment /i; = for 6 = 3, 

10, 30, 100, from (7) and (8). 








